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Abstract 

The X-ray background (XRB) is produced by a large number of faint sources distributed over a wide 

range of redshifts. The XRB carries information on the spatial distribution and evolution of these 

sources. The goals of the paper are: 1. to determine the redshift distribution of the soft X-ray background 

photons produced by all types of extragalactic sources, in order to relate fluctuations of the background 

to the large scale structures, 2. to determine the redshift distribution of the soft XRB produced by AGN 

OO ' in order to calculate the evolution of the AGN X-ray luminosity density. A set of major X-ray surveys is 

^ . used to determine the redshift distributions of the X-ray sources selected at various flux levels. Simple 

p^ I analytic fits to the data allow us to determine the smooth relationship between the redshift distribution 

and the source flux. The redshift distribution of the integral XRB flux is obtained by averaging the fits 

over the source counts. It is shown that the distribution of extragalactic XRB photons in the 0.5 - 2keV 

band is adequately represented by the function: dnxRB/dlogz. = 5.24 z.^'^^ exp(-z/0.63). The huge 

voids postulated to explain the cold spots in the CMB maps create dips in the total XRB flux. However, 

_H ■ the expected magnitude of the effect is comparable to the fluctuation amplitude of the XRB generated 

Q_|. by the individual sources contributing to the background. The cosmic evolution of the AGN X-ray lu- 



o 



s 



minosity density up to redshift of ~ 5 is calculated in an elegant and straightforward way. Systematic 



H ' uncertainties of the present method are assessed and shown to be small. At redshift greater than one the 

^ . present results could be compared directly with some recent estimates obtained in a standard way and 

the agreement between both methods is very good. 

KeywordsX-rays: diffuse background - intergalactic medium - X-rays: galaxies 

^ ■ 1 Introduction 
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The X-ray background (XRB) is generated mostly by discrete extragalactic sources (e.g. iLehmann et al. 



LCKgl 
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Q ■ 120011 : iKim et al.ll2007l. and references therein), predominantly by various types of Active Galactic Nuclei 



OO , (AGN) and cluster of galaxies. A question of flux, luminosity and redshift distributions of these sources 

^^ ' has been discussed in a great number of papers for the last 30 years. One of the major outcome of 

these investigations is th e conclusion that X-ray sources associa ted with the AGN are subject to strong 



S>^ . cosmic evolution (e.g. Miyaji et al.ll2000l. ISilverman et al.ll2007l. and references therein). As a result of 



the evolution, the redshift distribution of the XRB flux is wide. Thus, the integral XRB comprises the 
information on the large scale distribution of the X-ray sources over a wide redshift range. 

In the present paper the redshift distribution of the XRB photons is investigated in detail. The 
analysis is based on an extensive observational data selected from several published X-ray sky surveys. 
A convenient analytic approximations are applied to model the observed redshift histograms of the 
extragalactic X-ray sources selected at several flux levels. These distributions are weighted by the source 
counts and summed up to obtain the redshift distributions of the integral XRB. 

Next, the redshift distribution is used to define a relationship between the XRB signal and the large- 
scale fluctuations of the matter spatial distribution. The investigation has been rais ed by the recent report 
on the huge void generating a dip in the surface brightness of the radio background (|Rudnick et al.ll2007h . 



This void, responsible for the deficit of the radio surface brightness, allegedly generates also a cold spot 
in the CMB map via t he late -time integrated Sachs-Wolfe effect. Although the careful statistical analysis 



by ISmith & Hutererl (120081) has not confirmed the existence of this particular "cold spot" in the radio 
survey, a relationship between the large scale features of the matter distribution and the integrated sky 
brightness in various energy bands is a problem deserving some interest. 

Apart of the question of the XRB fluctuations induced by voids, the XRB redshift distribution is 
interesting per se, as it allows to assess the evolution of AGN phenomenon. A standard way to estimate 



a rate and type of this evolution is based on the examination of the X-ray luminosity functions deter- 
mined at the consecutive redshift bins. Unfortunately, the X-ray surveys produce flux-limited rather than 
luminosity-limited samples of sources. In effect, luminosity functions at diff'erent redshifts cover differ- 
ent luminosity ranges. This in turn severely impedes estimates of the luminosity function over a wide 
range of luminosities and redshifts. The total level of nuclear activity in galaxies within unit volume 
is given by the integral of the X-ray luminosity function. The question of the AGN cosmic evolution 
constitutes one of the central problems of observational cosm ology, and has been investigated for the 
last forty years (this issue was for the first time recognized by S chmidt II 1968 1). Here a question of the 
AGN evolution is addressed without the calculations of the X-ray luminosity function. The available 
observational data on X-ray source counts and redshifts are used to evaluate the redshift distribution as 
a function of source flux. This relationship and the source counts allow to calculate the redshift dis- 
tribution of the total XRB and the integral luminosity density generated by the AGN as a function of 
redshift. 

The organization of the paper is following. First, I present the formulae used in calculations. Next, 
in Sec|3l the basic information on the observational material extracted form the various archives is given. 
Since the comprehensive characteristics of the data and the source catalogs are described in the original 
papers, only the basic properties of the material are presented here. The numerical fits to the observed 
distributions are obtained in Sec. |4l In that section the calculations of the redshift distribution of the 
XRB photons are described in details. These results are applied in the Sec.|5]to quantify the relationship 
between the voids and the XRB variations. In Sec.|6]the distribution of the XRB flux produced by AGN 
is used to calculate the evolution of the AGN activity. Finally, potential sources of errors inherent in the 
present method are discussed in the Sec.|7J 

The 'canonical' standard cosmology is assumed throughout, with Hq = 70kms~'Mpc~\ D.^ = 0.3, 
and Qa = 0.70. 

2 Basic relationships 

In the present approach the X-ray source catalogs are used to construct redshift distributions of the 
extragalactic sources as a function of the source flux. At this stage, a question of the source (absolute) 
luminosities is not addressed. Let N{S) denotes the X-ray source counts, i. e. number of sources brighter 
than 5 in a unit solid angle, and fsiz) = dn{z.\S)/d log z is the redshift distribution of sources with flux 
S . Then, the redshift distribution of the XRB surface brightness, b{z), is equal to: 
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b(z)^- dSfsiz)S 



dN{S) 



dS 



(1) 



where the integration covers the entire "interesting" range of source fluxes S and b denotes the integral 
background flux: 

b= fb{z)dlogz= fdSS ^^^^■' 



dS 



(2) 



It is assumed that the fs (z) distributions are normalized: 

fs{z)dlogz = l. (3) 



!' 



Here the integration limits cover the total range of redshifts occupied by X-ray sources. The actual limits 
of the "interesting" range of fluxes is discussed below. 

The luminosity density, e(z), i. e. a total luminosity L generated in a unit comoving volume, V: 

dL 
s(z) = -. (4) 

is related to the flux distribution b(z) and the luminosity distance, Di{z): 

s(.z)-^=4nDl{z)b{z). (5) 

dlogz 



Table 1 : The X-ray surveys selected for the analysis 

Name Flux limits Number of sources 

(ergcm"^s"') All Extragalactic AGN 

RBS 1.0xl0-'2_5 0xl0-" 1764 1054 681 

NEP 5.0x10-14 -1.0x10-12 26\ 248 192 

RIXOS 2.5x10-14-5.0x10"" 393 318 235 

XMS 1.0x10-14-2.0x10-" 275 256 231 

CDFS 5.0x10-1^-1.0x10-" 205 201 197 

CDFN 1.5x10-1^-5.0x10-" 425 412 268 



The cosmological relatio nships betwee n the comoving volume and the luminosity distance in a flat 
space with A ?t is given by lHoggl (|l999h : 



dV c Dl 



dz Ho (1 + 2)2 ya^(i+z)3+Q^ 
Combining Eqs.|5]and[6]we finally get: 



(6) 



e(z) = 47r — — b{z) . (7) 

c InlOz 

Thus, to calculate the distributions b{z) and e(z), the source counts N{S) and the functions n{z\S) 
have to be determined using the observational material. In the next section the available X-ray surveys 
are examined from this point of view. 

3 Observational material 

Because the high imaging efficiency of X-ray telescopes in the soft band and numerous extensive iden- 
tifications programs, we concentrate on the X-ray band of 0.5 - 2.0 keV. At these energies a fraction 
of th e background resolved into discrete sources exceeds 90 % and is higher than in the other bands 



(e.g. iMoretti et al.ll2003L iBrandt & Hasinger II2005I) . Also a fraction of identified objects with measured 
spectroscopic or photometric redshifts is relatively high. 

Equation [T] shows that the b{z) distribution is sensitive to sources which perceptibly contribute to 
the XRB. Consequently, one needs to calculate the n{z \S) functions over a quite wide range of fluxes. 
To achieve this objective I have examined numerous X-ray surveys and selected several major source 
catalogs for further analysis. The overall characteristics of those catalogs are listed in Table [T] A 
common name of the survey/catalog is given in column 1. From the each catalog, sources for further 
processing have been extracted within fixed range of fluxes defined on column 2. The numbers of all 
sources, extragalactic and AGN with known redshifts are given in columns 3, 4, and 5, respectively. 

Statistical requirements which have to be satisfied by the source samples to properly determine the 
n(z \S) distribution are different than those for the luminosity function calculations. The individual sam- 
ple should contain sources from possible narrow range of fluxes, but the sample has not to be flux limited. 
The sample provides unbiased estimate of the redshift distribution as long as the process of identification 
and redshift measurements does not introduce spurious correlation between flux and redshift. 

3.1 The ROSAT Bright Survey (RBS) 

The identification program of the brightest sources detected in the ROSAT All-Sky Survey, known as 
ROSAT Bright Survey, resulted in a sample of 2072 sources with the total count rate above 0.2 s"! 



(|Schwope et all 12000) . More than 99.5 % of sources in the final catalogue is identified. The survey 



covers high galactic latitudes {\b\ > 30deg). After the removal of the Virgo clusters and Magellanic 
Clouds regions, the catalog contains 2012 sources. In the energy band of 0.5 - 2.0 keV, 1773 RBS 
sources generate flux between 1.0 x 10-" erg cm-^s-i (hereafter cgs) and 5.0 x lO-n cgs. As one might 



expect, only for a small fraction of the RBS sources the redshifts are undetermined and relatively large 
number of sources is associated with galactic sources, mostly late type stars and cataclysmic variables. 

Since the RBS sample covers rather wide range of fluxes, it is useful for our purposes to divide 
it into several subsamples with narrow flux limits and to estimate the n{z \ S ) distribution for the each 
set separately. We define the bright source sample, labeled RBS(b), which contains sources with 4.0 x 
10"^^ < S < 1.0 X 10"^' cgs. Of 365 RBS sources in this flux range, 129 is identified with galactic 
objects and for the other 15 sources redshifts have not been measured. Eventually, the sample contains 
221 extragalactic objects with known redshifts. Nearly half of the sample, viz. 97 sources are identified 
with clusters of galaxies and normal galaxies. The median fluxjlj in the extragalactic subsample S^ - 
5.5 X 10"'2 cgs. 

The flux limits of 2.0 X 10"'^ < S < 4.0 x 10"^^ cgs have been adopted for the medium flux sample, 
RBS(m). Within this flux range the RBS comprises of 669 objects including 272 galactic stars. Among 
the 397 extragalactic sources, 233 are identified with AGN with known redshifts. The median flux in 
the RBS(m) extragalactic sample, S^ - 2.7 x 10~'^ cgs. The faint source sample, RBS(f), contains 
sources with 1.0 x 10"^^ < S < 2.0 x 10~'^ cgs. Of 599 RBS sources in this flux range 268 is identified 
with galactic objects and 24 has unknown redshifts. Thus, the RBS(f) sample contains 307 sources 
extragalactic objects with the median flux of 1.6 x 10~^^ cgs. 

The redshift distributions of sources in the RBS (b), (m), and (f) samples are shown in three upper 
left panels in Fig. \T\ The integrals of the histograms for all the samples in Fig. [T] are normalized to 
unity. The distributions are plotted using logarithmic redshift bins with Alogz = 0.125. Each histogram 
is labeled with the survey name and the median flux in cgs. Analytic fits will be discussed in the next 
section. 

3.2 The ROSAT North Ecliptic Pole Survey (NEP) 

The deepest expo sure of the ROSAT All-Sky Survey (RASS) is centered at the north ecliptic pole 
(|Voges et al.lll999r) . The RASS of this region has be en used to construc t statistically well defined sample 



of X-ray sources above a flux limit ~2 x 1 '^ cgs (iHenry et al.ll2006h which have been foUowed-up by 



the optical observations (|Gioia et al.ll2003h . The identification rate in the final catalog of 443 sources is 
very high (99.6 %). 

Within the flux limits of 5.0 x 10"^^ and 1 .0 x 10"^^ cgs the NEP survey provided 361 sources. After 
excluding 113 galactic stars, we are left with 248 extragalactic sources; for 3 sources the redshift is 
unknown. The redshift histogram of 245 sources (including 53 clusters) is shown in Fig. [ij the median 
flux in this sample is equal to 1.2 x 10"'^ cgs. 

3.3 The ROSAT International X-ray/Optical Survey (RIXOS) 

This ROSAT medium-sensitivity survey consists of sources found in 82 PSPC pointing observations at 
high galactic latitudes (\b\ > 28 deg). A flux limit of 3 x 10~^^ cgs was adopted in 64 fields and 8 x lO"^'* 
cgs in the remainin g 18 fields. The so urce selection procedures, optical identifications and the final 



catalog are given by iMason et all (|2000l '). 

For the purpose of the present analysis, 393 sources with fluxes between 2.5 x 10"^"* and 5 x 10~^^ 
cgs have been selected. Within these flux limits 75 sources are associated with galactic stars. Of the 
remaining 318 sources, the redshifts of three objects are unknown, and 49 sources are still unidentified. 
The redshift distribution of 266 sources (including 33 clusters) is shown in Fig. [T] The median flux in 
this subsample S^ = 5.2 x lO"''* cgs. 

3.4 The XMM-Newton serendipitous survey (XMS) 

The XMM-Newton serendipitous survey (XMS) has been constructed in a similar way as the RIXOS. 
More than 300 sources have been isolated in 25 high galactic latitude {\b\ > 22 deg) pointings covering 



'The samples extracted from the RBS span over narrow ranges of fluxes and there is no major difference between the mean 
and median values. In some other samples investigated in this paper the median flux is distinctly smaller than the average. In 
those cases, the median flux is adopted as the argument in the n(z \S) functions. 
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Figure 1 : Histograms with labels - the redshift distributions of extragalactic sources in 9 samples con- 
structed in the present investigation. The data are binned with Alogz = 0.125; The integrals of all the 
histograms are normalized to unity. The sample designations and the median fluxes in cgs are given in 
the upper left corners; dotted curves - best 3 parameter fits for the each sample separately; solid curves 
- the best model for all the samples. The lower right panel: solid curve - the model redshift distribu- 
tion for the whole XRB (the normalization of the distribution is the same as for the histograms); dotted 
curve - analytic approximation of the model distribution; points - the distribution obtained by the linear 
interpolation between the raw histograms (see text for details). 



~ 3deg^ of the sky (IBarcons et al.ll2007h . In the 0.5 - 2.0 keV band the sample is complete above 
1.5 X lO"''* cgs, and contains many weaker sources. 

I have extracted from the original catalog 275 sources with fluxes between l.Ox 10"^^ and 2.0x 10"'^ 
cgs. The sample is completely identified; it contains 19 stars and 2 clusters of galaxies. For 23 objects 
the redshifts are unknown. The redshift distribution of 233 sources is shown in Fig.[T] The median flux 
in this subsample ^m = 1.9 x 10~'^ cgs. 



3.5 The Chandra Deep Field-South (CDFS) 



The 1 Ms Chandra observations known as the Chandra Deep Field South are described by lGiacconi et al. 
(|2002h . The catalog of sources detected in this field by two independent algorithms contains 304 objects, 
of which 275 have determined fluxes in the 0.5 - 2.0 keV band. For further processing 205 sources with 
fluxes in the range 5.0 x 10"'^ - 1.0 x 10 ~^^ cgs have been selected. Four sources ar e identified with 
stars. The redshif ts either spectroscopic (|Szokoli et all l2004i iRavikumar et all l2007h or photometric 
([Zheng et al.l l2004) are known for 200 sources; one source remains unidentified. 



The redshift distribution is shown in Fig.[T] The median flux in this sample S^n - 2.1 x 10"'^ cgs. 
The sample contains o ne galaxy group; tvy o other sources apparently are not associated with the activity 
in the galactic nuclei ([Lehmer et al.ll2006r) . It is likely, however, that more objects in the CDFS survey 
should be classified as ofF-nuclear sources. The available data do not allow for unambiguous separation 
of AGN and off-nuclear sources at the low flux levels in the CDFS. This reservation holds also for the 
CDFN samples below. 



3.6 The Chandra Deep Field-North (CDFN) 

The ultra deep Ch andra field, 2 Ms exp osure, CDFN, resulted in a catalog of 503 sour c es de tected 
over 0.12 sq. deg. ([Alexander et al.ll2003h . Optical follow-up observations by iBarger et al.l (l2003h have 
rend ered a large number of spectroscopic and ph otometric redshifts. Several more redshifts are taken 
from iReddv et all (l2006h . bonlev et al.r(l2007h and ICeorgakakis et all (l2007h . 

In the present investigation, the CDFN catalog has been divided into two samples of bright (b) 
and faint (f) sources. The (b) sample contains 181 sources between 2.5 x 10~'^ and 5.0 x 10~'^ cgs; ten 
sources have bee n identified with the galactic stars; at least one has been categorized as 'starburst' galaxy 
(IGeorgakakis et al. 2007 ). for 49 objects the redshifts have not been measured. The sample consists of 



122 sources, mostly AGN. The median flux in this subsample Sm = 7.3 x 10^ cgs. 

The faint CDFN sample has been selected between 1.5 x 10"^^ and 2.5 x 10~'^ cgs. Among 244 
sources satisfying th ese flux limits, three sources have been identified with stars, 5 1 - with the starburst 
galaxies (Georgakakis et al.ll2007) and for 43 objects the redshifts have not been measured. The final 
sample used in the calculations contains 198 sources with 147 confirmed AGN. The median flux in the 
sample S^n = 9.0 x 10"^^ cgs. The redshift histograms for the (b) and (f) samples are shown in Fig.[T] 

The numbers of objects unidentified or without redshift are in some samples quite large. Hence one 
could expect that the corresponding redshift histograms are not representative for the whole population 
of sources at given flux. Below this question is de facto worked out where we construct an analytic 
function which simultaneously fits all the histograms. 



4 Approximations and fits 

The redshift distribution of sources selected at fixed flux, n{z\S), is a intricate function of a number of 
parameters, such as the luminosity function, the relationship between the luminosity and observable flux, 
and the relationship between the volume and redshift. The luminosity function itself depends on redshift 
and both the latter relationships depend on the cosmological model. However, the existing estimators 
of the n{z\S ) function represented by the nine histograms in Fig. [U are strongly affected/degraded by 
the statistical noise. It implies that a simple analytic function with 2-3 free parameters will provide a 
statistically satisfactory fit to the observed distributions. 

It appears that the histograms in Fig. [T] are adequately reproduced by: 



fsiz) = fofe- 



(8) 



0.1 



1 III 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 


1 1 i llli 1 i 1 1 1 1 III 1 II 


1 1 1 1 1 1 1 1 1 1 II 


i 


i 




- 


i 




5, 






\ 




' 




X 


" 


X X 


X 


X 
X 


X X 
X 


" 






'i. 


; 






■■i 


- 






\ 











10 



10-17 10-is 10-'^ 10-'* 10-13 



10- 



10- 



Figure 2: Distributions of the best fitted parameters a and Zc vs. the median flux, S^- Points with error 
bars - Zc (left-hand ordinate), crosses - a (right-hand ordinate). The assumed Zc ~ S relationships used 
to model the XRB redshift distribution is shown with the dotted lines (see text for details). 



where /o = foiS), a = a{S), and Zc = ZciS) are three parameters fitted to the the histograms n{z\Si), 
i = 1, ..., 9. In Fig.[T]the least square fits for all the distributions are shown with the dotted curves. Apart 
from a few pronounced features visible in the plot s which represent the large scale structures reported 
in the literature (e.g. Barger e t al.ii2002 . iGilli et al.ll2003i) . analytic fits seem to adequately reproduce the 
observed distributions. 

It is found that only Zc is strongly correlated with S^, while the fits do not indicate any statistically 
significant correlation between a and S^- In Fig.[2]the best fit values of a are shown with crosses. The 
labels and scale on the left-hand ordinate refer to Zc, and on the right-hand - to a. Since the simultaneous 
fitting of a and Zc introduces a spurious correlation between these two parameters, the a parameter has 
been fixed at the average value found for the 9 samples, a = 1.934. Effectively, it means that the shape 
of the n{z\S ) function is fixed and the only dependence on S is limited to the horizontal shift along the z 
axis. In Fig. |2]the best fit parameters Zc found for the fixed a are plotted against the median flux ^m- In 
agreement with the expectations, the Zc increases with diminishing flux S^ over a wide range of fluxes. 
However, a clear flattening of the relationship is observed below ~ 10"'^ cgs. This apparent absence of 
correlation between Zc and S results from the well-known fact that in the X-ray surveys at low flux levels 
the maximum detected redshift remains stable while significantly increases fraction of intrinsically weak 
sources. 

The points in the Zc - S^ relationship above ~ 10"^^ cgs seem to follow the power law. In the 
subsequent calculations it is assumed that this relationship is in fact well approximated by the power law 
in the whole range of fluxes between 10~^^ and 10~^^ cgs, although the data coverage is rather sparse. 
Below S w lO"^'^ cgs the data are insufficient to delineate precisely the Zc - S^ relationship. I have 
assumed tentatively that Zc remains constant and is equal to the average value found for three Chandra 
samples. The Zc in CDFN(b), CDFS and CDFN(f) are equal to 0.69, 0.69 and 0.51, respectively. The 
average weighted by the uncertainties Zc = 0.61. In Fig. |2]the model relationship Zc - S used in the 
calculations is shown with the dotted line. 

The Zc ~ S relationship with the fixed a parameter and fixed normalization of the integral: 

Jfs{z)dlogz = l, (9) 

eliminates formally any free parameters in the fitting the analytic function to the nine histograms in 



Fig. 121 

fs(z) = ^z'^e-^'-, (10) 

where a = 1.934, T(a) = 0.9739 is the gamma function and z.c is specified for each sample by 
the Zc ~ S relationship. Analytic distributions defined in Eq. [TO] are shown in Fig. \T\ with the solid 
curves. In all the histograms the model distribution is astonishingly close to the corresponding best 
three-parameter fit represented by the dotted curves. Most deviations visible in some plots are easily 
explained by the statistical natur e of the problem and/or the large sca le structures present in the catalogs 
based on the localized sky area ("Bar ger et al.ll2002L iGilli et al.ll2003b . Systematic shifts between the fits 



are present in three histograms below S = 10"^"* cgs. It is a direct result of the assumption of a single 
Zc - 0.61 value for all three Chandra samples. It is noticeable that the constant width (in logz) model 
fits adequately represent the data over the full range of fluxes. Small differences in the width between the 
three-parameter fits and the final model which are visible in the NEP, RIXOS and XMS data, apparently 
do not represent the systematic effects. In the NEP and RIXOS histograms the final model is slightly 
narrower than the individual fits, while in the XMS sample it is wider. 

To effectively use the Eq. [T]one needs the represen tation of the s ource counts dN(S)/dS over the 



whole range of fluxes S . The parametrization by iMoretti et all (l2003h adequately suits the present cal 



culations. The smooth functional form for A'^(5' ) proposed by Moretti et al. accurately reproduces the 
observed counts below 10~'^ cgs down to Chandra threshold of ~ 2 x 10^^^ cgs. Sources within these 



flux limits generate more than 90 % of the XRB and smooth extrapolation of the IMoretti et all (|2003h 
counts down to ~3 x 10"'^ cgs is consistent with the entire XRB. Substituting all the components into 
Eq. [T]we finally get the redshift distribution of the XRB photons. It is shown with the solid curve in the 
bottom right panel in Fig. \T\ The same normalization has been applied to facilitate comparison with the 
distributions derived for the individual samples. Points in the plot are discussed below in the Sec.|7j 

A suitable representation of the b{z) distribution has been found using a smooth function of the same 
form as for the individual redshift histograms. The function: 

Mz) - 5.24 zi-52 ^--'/o.63 ^ (11) 

reproduces the derived distribution of b{z) with the relative error of less than 4 % for 0.06 < z < 6. It 
is shown with dotts in the bottom right panel (with normalization rescaled to conform to all the plots in 
Fig.©. 

5 XRB and Supervoids 

The distribution of the XRB photons b(z) peaks at redshift z * 1 and 50 % of the background originates 
between the redshifts of 0.4 and 1.4 (for 80 % the redshift limits are 0.2 and 2.1). Thus, very large struc- 
tures of the matter distribution at redshift within these limits would generate fluctuations of the integral 



XRB. As an example I discuss below the X-ray signature of the huge void postulated by iRudnick et al. 



(120071) . The arguments based o n the radio survey in fa vor of the void with a radius of ~ 140 Mpc in 



Eridanus have been questioned r Smith & Huterenl2008h . Nevertheless, the Integrated Sachs-Wolfe ef 



feet operating on extremely large structures of matter remains a valid explanation of the strongest CMB 
fluctuations. 

The redshift separation, Az, correspo nding to the far side and the near side of the void with diameter 
Rg centered at redshift z is equal to (e.g. iHogg II 19991) : 

op 

Az = —^ V^m(l + z)' + Qa ■ (12) 

Using the b(z) distribution we assess the fractional deficit of the XRB 6 = Ab/b created by the completely 
empty region of size 280 Mpc. Such void would generate \6\ = 5.7 % and 6 = 4.9 % at redshifts z = 0.5 
and 1, respectively. Assuming spherical shape of the void, its angular diameter would be 8?4 and 4?7 
at these redshifts. The XRB depression produced by the void should be compared to the XRB intrinsic 
fluctuations resulting from the discrete nature of sources generating the background. Assuming purely 



random distribution of sources, the rms fluctuations of the XRB are defined by the source counts N{S): 

1/2 



o-b 



f 



dS S 0) 



dNiS) 



dS 



(13) 



where co is the solid angle subtended by the investigated area. Using the the lMoretti et all (120031) counts 
and 5 max - 1 • 10"'' cgs (the amplitude of the XRB fluctuations is dominated by the contribution of 
sources at the bright end of counts), we get 6b/b = 0.035 and 0.024 for the circular areas of radius 2° 
and 3°, respectively. Thus, at z = 0.5 the signal-to-noise ratio for the void detection amounts to ~2.4. In 
the case of z = 1 the S/N drops to just 1.4. To reduce the ampUtude of the XRB fluctuations one should 
remove from the XRB the contribution of bright sources. If the S^ax is decreased to 1 • 10~'^ cgs, the 
significance of the void signal reaches 3.5 cr at z = 0.5 and 2.0 cr at z = 1. So, only the low redshift voids 
would produce the XRB deficits significantly stronger than the statistical fluctuations. 

6 Redshift distribution of the XRB and the AGN evolution 

To assess the distribution of the XRB generated just by the AGN, I have repeated all the procedures 
described in the previous sections using the samples constructed exclusively from the AGN. The AGN 
sources are easily separated from the clusters and nearby normal galaxies. However, the distinction 
between the nuclear activity and stellar emitters in the case of distant and weak sources becomes prob- 
lematic. Such sources are present in both Chandra surveys. One should keep in mind this limitations 
in the present investigation. Nevertheless, the AGN are a dominating constituent of all the samples ex- 
ploited in the paper and even the moderate contamination of the AGN subsamples with the off-nuclear 
sources would not affect significantly our calculations (see below). The numbers in the AGN samples 
are smaller than in the full samples and parameter estimates are subject to slightly larger uncertainties. 
Clusters and normal galaxies populate on the average lower redshift bins and the histograms for the AGN 
analogous to those in Fig. [J are shifted towards the higher redshifts. We notice also a weak correlation of 
the best fit a parameter with the source flux - a shape of of the redshift distribution (in log z bins) varies 
with Sm- In Fig. [3] the values of a in the nine samples are shown with crosses. The regression line of 
log a on log S is used to fix the value of a for the each sample and to calculate the best fit parameter Zc- 
These new Zc are shown in Fig. [3] with the squares. Finally, the best fit line logZc ~ log ^m is calculated 
for six brighter samples to obtain Zc for 5 > 1.9 • 10~'^ cgs. A constant Zc is assumed for lower fluxes, 
and the complete Zc ~ S^ relationship adopted for further computations is shown in Fig. [3] with the 
dotted lines. 

The X-ray source counts used in the present case should be limited to the AGN only. Two other 
major classes of sources contributing to th e counts are associa ted with clusters and normal/starburst 



galaxies. The analytic formula obtained by 'Moret ti et al.l (120031 ) quite accurately represents counts of 



all the types of extragalactic sources, but the relative contribution of the each class in the total counts is 
not well established. One should notice, however, that most of the XRB is produced by sources in the 
middle range of fluxes considered here, while the cluster contribution is significant only at the bright end 
of counts and the normal and starburst galaxies populate mostly the faint end of counts. Vertical bars at 
the bottom of Fig.[3]divide the flux range of 10"'^ - 10"'' cgs into 9 contiguous bands corresponding 
approximately to fluxes surveyed by the source samples defined in the paper. The numbers between the 
bars give the relative contribution of each flux band to the total XRB. 

To extract the cluster and starburst galaxies contributions we coiTccted the total counts in the fol- 
lowing way. In the RBS(b) sample clusters constitute 40% of all the extragalactic sources. The slope 
of the cluster counts at the bright end amounts approximately to 1.3 (|De Grandi et all 1 19991) . Substan- 



tially flatter slope than that for the total counts reduces the relative cluster contribution at lower fluxes. 
Although the cluster counts are not well constrained below ~ 10"'^ cgs, their contribution to the total 
counts drops at the faint end of counts to a negligible level. The source counts a ttributed to th e AGN 



are assessed by subtracting the cluster counts from the total counts defined by the IMoretti et all (|2003h 



formula. The normal and starburst galaxies are relatively abundant in the CDFN(f) sample. Of 241 ex- 
tragalactic sources, 147 have been classified as 'AGN', 51 as 'starburst' and for the other 43 the redshift 
is unknown. The absolute maximum content of the non-AGN sources in the CDFN(f) sample amounts 
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Figure 3: Distributions of best fitted parameters a and Zc vs. the median flux, S^, in the AGN sam- 
ples. Points with error bars - Zc (left-hand ordinate), crosses - a (right-hand ordinate). The dashed line 
represents the regression line for the a ~ Sj-n relationship; the Zc ~ S function used to model the AGN 
contribution to the XRB is shown with the dotted lines (see text for details). The vertical bars at the 
bottom mark the effective flux limits of the analyzed samples; numbers indicate a relative contribution to 
the XRB generated by the AGN within these limits. 



to (51 + 43)/241 = 39%, assuming that all sources with undetermined redshift are starburst galaxies. 
The amount of the non-AGN sources at higher flux levels drops quickly. In the XMS sample none ex- 
tragalactic source with known redshift has been classified as normal or starburst galaxy. The maximum 
possible contrib ution of the starburst galaxies at the low flux levels has been accounted for by flattening 
the slope of the Moretti et al.l (120031) counts below 10"^"* cgs to reproduce the reduction of the AGN in 
the CDFN(f) sample by 39 %. The counts modified this way have been substituted into Eq. [T]to obtain 
the redshift distribution, Z7agn(z)- 

One can use the observed distribution of the background flux produced by the AGN to calculate the 
cosmic history of the luminosity density generated by these objects. The ^agn(z) distribution is inserted 
into Eq.|7] which relates the cosmological evolution of the X-ray luminosity density, e(z), to the redshift 
distribution of the background, b(z). Variations of the luminosity density obtained this way are shown in 
Fig. |4] with the solid curves. The data are displayed in three panels as a function of redshift, z, logarithm 
of (1 -I- z) and the cosmic time, assuming fo = 13.47 • 10^ years for the present age of the Universen. 

The accuracy of the present e(z) estimate depends strongly on a quality of our b{z) fits. Relatively 
small numbers of sources at redshifts below ~ 0.03 and above ~ 3 generate large statistical fluctuations 
and weakly constrains the analytic fits b{z) in these redshift ranges. Hence, the present estimates of e{z) 
are also subject to large uncertainties at low and high redshifts. 

In order to assess the importance of the A'^(5') uncertainties on the pr esent estimates of e( z), I have 
plotted in Fig. |7] with the dotted curves the s(z) function using the original iMoretti et all (120031 ) formula, 
i.e. assuming no corrections for clusters and starburst galaxies. The discrepancies between both solutions 
do not exceed 20 % for redshifts below ~ 3. It implies that our procedure to isolate the contribution of 
AGN from the total counts, albeit crude, does not contribute significantly to the final errors of e(z). 



^For the cosmological model defined in Sec.[T]and using the formulae given bv lHoggl (119991) 
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Figure 4: The cosmic evolution of the X-ray luminosity density generated by AGN as a function of: (a) 
- redshift z, (b) - log(l +z) and (c) - the cosmic time. The plots are based on the soft X-ray emission only 
and do not include radiatio n absorbed and r e-emi tted in different energies. Points with the error bars in 
the panel (a) are taken from lHasinger et al.l (l2005h and represents only the type-1 AGN. 
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7 Discussion 

The main objectives of the present investigation, viz. estimates of the redshift distribution of the XRB 
photons, b{z), and the evolution of the AGN luminosity density, e(z), have been achieved using the 
smooth, analytic fits to the observed source redshift histograms. The present method is conceptionally 
simple and computationally straightforward. Unfortunately, it does not provide error estimates. The 
major sources of uncertainties have been indicated in the previous section. Here a quantitative estimate 
of the errors is discussed. 

The errors of the present measurement of b{z) are generated by the statistical nature of the inves- 
tigated material and a chain of approximations applied to substitute the observed redshift distributions 
centered on a selected fluxes by an analytic function fs (z) continuous in both parameters, z and S . In fact, 
the visual inspection of the analytic fits displayed with the solid curves in Fig.[T]reveals some deviations 
from the redshift histograms. To estimate the significance of these differences, the calculations have 
been performed using the actual histograms shown in Fig. [T] with broken solid lines instead of fs(z)- For 
the each value of flux S in the range 10^'^ - 10"^' cgs, the corresponding redshift distribution has been 
obtained by the linear interpolation between two histograms from the samples centered on the median 
fluxes nearest to 5 . The results of this procedure are shown in the lower right panel of Fig. [T] with dots. 
Generally good agreement between the distribution of points and the solid curve proves that the analytic 
approximations do not introduce perceptible systematic errors in the present investigation. It appears 
that the relatively large deviations for three data points (centered at redshifts: 0.087, 0.65, and 1.16) 
result purely from the large scale structures. This is particularly likely for the first bin (0.075 < z < 0.1), 
where the discrepancy between the fits is produced entirely by the excess of sources in the localized NEP 
survey. 

One should also notice, that the uncertainties of our main results are only weakly affected by the 
limited statistics of the individual samples and histograms. This is because the final distributions are 
obtained by averaging the individual distributions and this procedure effectively reduces statistical fluc- 
tuations. 

In the upper panel of Fig. |4]the AGN emissivity calculated by 'Hasinger et al. ( 2005h is shown. The 



points with the error bars are redrawn here from their original paper. Hasinger et al.l ( 12005) apply more 



stringent criteria to select sources and use only well defined samples of type-1 AGN. Optically these 
objects are identified by the broad Balmer emission lines, while using the X-ray criteria, they have 
unabsorbed spectra indicating low intrinsic column densities. In the present analysis I have included all 
objects in which th e X-ray emission origi nates in the active nuclei. Thus, our results cannot be compared 



directly to those by iHasinger et all (120051) . Nevertheless, despite entirely different method applied in the 
present paper, the distributions show good agreement at redshifts above ~ 1. Although most of the 
apparent discrepancies, which reach a factor of 3 at z « 0.5, are probably due to the distinct selection 
criteria of both investigations, one cannot exclude that some differences are caused by unrecognized 
systematic effects inherent in one or both methods. 

The present method of the luminosity density calculations has also some disadvantages. In our 
approach the absolute luminosities of the individual objects are not determined. Consequently, only the 
integral luminosity density is obtained, and the cosmic evolution of any selected AGN luminosity class 
has to be studied by means of the standard methods. 
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